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ABSTRACT 

Using archival Hubble Space Telescope (HST) imaging data, we report the multiband photometric 
properties of 13 ultraluminous X-ray sources (ULXs) that have a unique compact optical counterpart. 
Both magnitude and color variation are detected at time scales of days to years. The optical color, 
variability, and X-ray to optical flux ratio indicate that the optical emission of most ULXs is dominated 
by X-ray reprocessing on the disk, similar to that of low mass X-ray binaries. For most sources, the 
optical spectrum is a power-law, F v oc v a with a in the range 1.0 to 2.0 and the optically emitting 
region has a size on the order of 10 12 cm. Exceptions are NGC 2403 X-l and M83 IXO 82, which 
show optical spectra consistent with direct emission from a standard thin disk, M101 ULX-1 and M81 
ULS1, which have X-ray to optical flux ratios more similar to high-mass X-ray binaries, and IC 342 
X-l, in which the optical light may be dominated by the companion star. Inconsistent extinction 
between the optical counterpart of NGC 5204 X-l and the nearby optical nebulae suggests that they 
may be unrelated. 

Subject headings: black hole physics — accretion, accretion disks — galaxies: stellar content 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are nonnuclear 
X-ray sources with luminosity above the Eddington limit 
of stellar mass black holes (Lx > 3 x 10 39 erg s _1 ) as- 
suming isotropic emission. ULXs with fast variability 
or long-term chaotic flux variation may contain accret- 
ing black holes in an attached binary system. They 
could be powered by accretion onto either intermedi- 
ate mass black holes of 10 2 — 10 4 M Pl at sub-Eddington 
rate (iColbert fc Mushotzkvlll999l ; lMakishima et~aT1l2000l: 
IKaaret et al.l I200lh or normal or slightly massive stellar 
black holes of < 100 M Pl at super or n ear Eddington rate 
(|Watarai et aUl2001b IBegelman|[2002l ). 

In black hole binaries, the X-ray emission probes the 
deep potential at regions close to the event horizon 
while the optical emission is important to constrain their 
masses, interaction with environment, and evolutionary 
history of the binary. Optical observations of ULXs have 
played a key role in unveiling the nature of these sources. 
Many bright ULXs are found to be spati ally associated 
with an optical nebula (jPakull et al.ll2003l ). These nebu- 
lae are large with an extent of tens to hundreds pc, and 
expanding at a velocity in the order of 10 2 km s _1 . Most 
of them are powered by strong shocks created by out- 
flows from the binary hitting the surrounding medium, 
and in a few cases, extraordinary photoionization is re- 
quired (jKaaret et alj|2004 IKaaret fc Corbel l2009). The 
photoionized nebulae can be used as bolometers to ex- 
treme ultraviolet and soft X-ray photons, and suggest 
that the ULX luminosities are indeed high. Point-like op- 
tical objects as compact counterparts of ULXs have been 
identified for 10 sources in the literature. The companion 
star, intrinsic emission from the disk at large radii, and 
disk irradiation may all contribute to the optical emission 
and even dominate at different wavelengths. Binary pop- 
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ulation synthesis indicates that colors and magnitudes of 
the optical counterpart of ULXs can help distinguish be- 
tween stellar mass and intermediate mass black holes, 
and a sample of 6 sources f avors the latter interpretation 
(Madhusudhan ct al. 2008). The continuum component 
in the optical spectrum of NGC 5408 X-l is consistent 
with that from a standa rd accretion disk with irradiation 
(|Kaaret fc Corb el 2009). However, this result is obtained 
from a ground-based telescope so the nebular contribu- 
tion cannot be directly resolved and has to be subtracted 
based on model prediction. Therefore, in most cases, 
HST observations are required for the study of the com- 
pact optical counterpart of ULXs, in particular for the 
continuum emission. 

There are 10 ULXs with a unique optical counterpart 
identified in the literature. We make 3 new identifica- 
tions in this paper. Here, we analyze all available HST 
imaging data for these 13 ULXs, see Tabled to investi- 
gate their multiband spectra and variation at long time 
scales. The observations and data analysis are described 
in §[21 and the results are discussed in §[3] and conclusions 
summarized in § 2] 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Sample 

Up to beginning of analysis work for this paper, there 
have been 10 ULXs with a uniq ue optical counter - 
part identified: Hol mberg II X-l (IKaaret et al.ll200l . 
Holmberg IX X-l (IGrise et all 120111), IC 342 X-l 



dFeng fc Kaaretll2008l) . M81 ULS1 (ILiu et alJl200l). M8 1 
X-6 (ILiu et al.ll2002D. M101 ULX-ilKuntz et al.112 005) , 
NGC 1313 X-2 (|Ramsev et all 120061: IPakull et al 2009: 
ILiu et alJl2007l). NGC 5204 X -l (ILiu et all 120041) . NGC 
5408 X-l (ILang et alJ l200l. a nd NGC 6 9 46 U LX-1 
(IKaaret et al.M2010D Moreover. IPtak et al.l (|2006D and 
iRoberts et al.1 (|200c? ) reported possible optical counter- 
parts of another 50 ULXs. However, these sources ap- 
pear to have more than one optical counterparts shown 
in the X-ray error circle. Here, using Chandra and HST 
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TABLE 1 

Positions of the reference objects, ULXs, and their 
unique optical counterparts 



TABLE 2 

List of the ULXs studied in this paper. 



source 


instrument R.A. 


decl. 


A 






(J2000.0) 


(J2000.0) 


(") 


NGC 4559 X-7 


Ref. 1 


Chandra 


12 36 00.353 


+27 54 56.61 


0.271 




HST 


12 36 00.354 


+27 54 56.67 




ULX 


Chandra 


12 35 51.724 


+27 56 04.35 


0.052 




corrected 


12 35 51.725 


+27 56 04.41 


0.28 


counterpart 


HST 


12 35 51.719 


+27 56 04.41 




M83 IXO 82 


Ref. 1 


Chandra 


13 37 14.764 


-29 54 28.47 


0.318 




HST 


13 37 14.761 


-29 54 28.26 




Ref. 2 


Chandra 


13 37 19.986 


-29 53 18.09 


0.169 




HST 


13 37 19.979 


-29 53 17.68 




ULX 


Chandra 


13 37 19.804 


-29 53 48.80 


0.029 




corrected 


13 37 19.799 


-29 53 48.43 


0.18 


counterpart 


HST 


13 37 19.797 


-29 53 48.58 




NGC 2403 X-l 


ULX 


Chandra 


07 36 25.567 


+65 35 39.87 


0.6 




Chandra 


07 36 25.558 


+65 35 39.98 


0.6 




mean 


07 36 25.563 


+65 35 39.93 


0.42 


counterpart 


HST 


07 36 25.562 


+65 35 40.08 





Note. — The used Chandra observation is 2686 for NGC 4559, 
793 for M83, and 2014/4630 for NGC 2403. The used HST datasets 
are j92w02030, j9h811010, and j9h803010, respectively for the three 
targets. A is the 90% statistical error radius; for NGC 2403 X-l, it is 
the absolute astromctry error of Chandra. 



ULX 


D 






E n 


Neb. 


Ref. 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


Holmbcrg II X-l 


o.UO 


U.Uoz 


U 


U 1 


v 
1 


1 A 
1-4 


Holmberg IX X-l 


3.6 


0.079 





26 ± 0.04 


Y 


1,2,5,6 


IC 342 X-l 


3.3 


0.565 





82 


Y 


1,2,7,8 


M81 ULS1 


3.63 


0.080 






N 


1,5 


M81 X-6 


3.63 


0.080 






Y 


1,2,5 


M83 IXO 82 


4.7 


0.066 






N 


1,9 


M101 ULX-1 


7.2 


0.009 





13 


Y 


1,10,11 


NGC 1313 X-2 


3.7 


0.085 





13 ± 0.03 


Y 


1,2,9,12 


NGC 2403 X-l 


3.2 


0.040 






N 


1,13 


NGC 4559 X-7 


10 


0.018 






N 


1,9 


NGC 5204 X-l 


4.3 


0.013 






Y 


1,2,14 


NGC 5408 X-l 


4.8 


0.069 





08 ± 0.03 


Y 


1,15-17 


NGC 6946 ULX-1 


5.1 


0.343 





co+0.08 
ou -0.07 


Y 


1,18-22 



Note. — D is the adopted distance to the host galaxy in 
Mpc. Eg is the Galactic extinction and E n is the total extinction 
estimated from the surrounding optical nebulae. Column (5) in- 
dicates whether or not the source is spatially associated with an 
optical nebula. Column (6) lists the referen c es from where the 
i nformation is obtai ne d: fechlegel et al.l fl99l; ^ Pakull ct aQ 120021 ; 
Stewart et al.l 120001 ; IHoessel et alj 119981; feeedman et al.1 11994 
Grise et al.l I201U; IGrise et al.l 120061: * ISaha et al.l 120021; ITullv 



19.S.S 



2008; 



Abolmasov et al.l (200^ 1J |Stetson et al.l I1998t IGrise et al 
fieedman fc Madorel [1988; ^ITullv. Shava fc Pierce! 1199 ~ 
E Karachcntscv ct al. 2002; "Kaarct & Corbel 2009; "Pakull ct al 
I2003t it lde Vaucouleurs. G.I 1197a ^lAbolmasov et alj [200 
^Blair et al.l 1200 ll ^ IRoberts &; Colberti 120031 : ^ IBlair fe Fese 



11994 



data up to date, we surveyed ULXs in these catalogs for 
possibilities of improving the relative astrometry, which 
leads to the identification of a unique counterpart for an- 
other 3 sources, M83 IXO 82, NGC 2403 X-l, and NGC 
4559 X-7. 

The details about the astrometry correction have been 
described in iFeng fc Kaaretl ([200 8). The relative as- 
trometry between Chandra and HST is improved by 
aligning objects (hereafter referred as reference objects) 
shown on both. For Chandra, deepest observations 
with the Advanced CCD Imaging Spectrometer are used. 
The X-ray positions are calculated using the wavdetect 
tool in CIAO from exposure-corrected images in the 
energy range of 0.3-8 keV. The HST images are first 
registered to th e Two-Micron All Sky Survey (2MASS; 
ISkrutskie et"all[2006l) source frame to achieve better ab- 
solute astrometry before searching for reference objects. 
For NGC 4559 X-7 and M83 IXO 82, there is 1 and 2 
reference objects, respectively, identified for image align- 
ment. The corrected X-ray error circle of the ULX af- 
ter alignment, taking into account of all uncertainties, is 
around 0.2"-0.3", leading to a unique identification (Ta- 
ble [J). For NGC 2403 X-l, although we cannot find a 
reference object to align the images, there is only one 
HST object shown around the mean position from two 
Chandra observations. We thus regard this object as the 
optical counterpart of the X-ray source. The HST images 
around the 3 ULXs are shown in Figure [T] The solid cir- 
cle indicates the 90% absolute position error of the ULX 
from Chandra with a radius of 0.6", and the dashed cir- 
cle indicates the aspect-corrected X-ray position relative 
to optical images with a radius listed in Table [TJ The 
uncertainty in the ULX error circle is calculated via er- 
ror propagation using the uncertainties in the position of 
the ULX relative to the reference source(s) in both the 



X-ray and optical. This alignment technique, based on 
a small number of reference sources, has been used for 
most ULX optical identifications. In several cases, the 
optical counterpart identification has been confirmed via 
the detection of high excitation emission lines that could 
not be produced by a star with the given optical col- 
ors, but are easily produced by X-ray illumination of an 
accretion disk. 

As a whole, a list of the 13 ULXs in our sample with 
related information, like the distance to their host galax- 
ies, Galactic extinction along the line of sight, their en- 
vironmental associations, and the references, is shown in 
Tabled! 

2.2. Photometry 

The HST data used in the paper are observations with 
either the Advanced Camera for Surveys (ACS) instru- 
ment or the Wide-Field Planetary Camera 2 (WFPC2). 
For ACS data, the photometry is performed on drizzled 
images produced from standard pipeline calibration us- 
ing the APPHOT package in IRAF with an aperture ra- 
dius of 0.15"-0.2". For Holmberg IX X-l, IC 342 X-l, and 
NGC 4559 X-7, each has a weak optical source located 
at about 0.2" from the counterpart of the ULX. To mini- 
mize the contamination, we used a smaller aperture with 
radius of 0.1"-0.15" for photometry and added the aper- 
ture correction error into the results. For WFPC2 data, 
the photometry is performed on the cOf images using 
HSTPHOT 1.1, which outputs aperture-corrected VEGA 
magnitudes. As some sources are shown in crowded re- 
gions, the choice of background may slightly alter the de- 
rived net flux, which is particularly prominent for sources 
associated with an optical nebula that a flux difference up 
to 30% could arise whether or not the nebular emission 
is subtracted (e.g. [Kaaret et al.l 120101) . Therefore, the 
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Fig. 1.— Upper panels: HST images around NGC 4559 X-7, M83 IXO 82 and NGC 2403 X-l . The solid circle indicates the original 
X-ray position of the ULX, with a radius of 0.6" corresponding to the absolute astrometry error of Chandra. The arrow points north and 
has a length of 1". For NGC 4559 X-7 and M83 IXO 82, the dashed circle is the corrected X-ray position after alignment. For NGC 2403 
X-l, the dashed circle is the mean position of the two Chandra observations. A unique optical counterpart is identified for each source. 
Lower panels: SEDs of the ULXs and nearby sources. The nearby sources are marked on the images in the upper panel. In all three 
cases, the SED of the ULX counterpart is consistent with a single power-law, while the nearby stars show a turn over at short wavelengths. 



uncertainty introduced by sky background subtraction is 
taken into account, estimated using various background 
regions or background estimate means. 

The total extinction along the line of sight to the source 
is adopted in two ways. For those associated with an op- 
tical nebula and the extinction to the nebula has been 
measured from emission lines, the nebular extinction is 
adopted (see Tabled- Otherw ise, we adopt Galactic ex- 
tinction (jSchleeel et al.l [19981) for reddening correction. 
The extinction can also be estimated from the neutral hy- 
drogen column density, which is derived via X-ray spec- 
tral fitting. For Holmberg II X-l, IC 342 X-l, M101 
ULX-1, NGC 1313 X-2, NGC 5408 X-l, and NGC 6946 
ULX-1, which reside in nebulae, the ratio of Galactic to 
nebular extinctions (E g /E n ) are 0.46, 0.69, 0.07, 0.65, 
0.86, and 0.69, respectively, and the ratio of X-ray to 
nebular extinctions (E^/E n ) are 1.0, 1.42, 0.85, 3.0, 1.63, 
and 1.42, respectively. For Holmberg IX X-l, the nebu- 
lar extinction is found to vary with position within the 
nebula from 0.09 to 0.26 (|Grise et al.l 1201 ID . and con- 
sequently causing E s /E n to vary between 0.3 and 0.9, 
and E x /E n between 1.3 and 3.9. It seems that E x of- 
ten overestimate the total extinction. Plus, the X-ray 
column density usually varies dramatically, causing ad- 
ditional difficulty in deriving the extinction. Thus, we 
adopt Galactic extinction when the nebular extinction is 
unavailable, and this should be treated as the lower limit. 
A Ga lactic extinction law (Ca rdelli. Clavton fc Mathisl 
1989) is assumed for reddening correction. Except for 
the small handful of images in the near UV (F330W fil- 
ter), the resultant magnitude varies very little, by 0.002 



typically, if we choose a different reddening law for the 
excessive extinction above the Milky Way, like that for 
the Small Magellanic Clouds where the metallicity is sig- 
nificantly lower. We hence use the Galactic reddening 
law for all objects despite that metal poor environments 
have been sugges ted for a few of them (e.g. lMirioni|[2002t : 
ISoria et all 12001 . 

A power-law spectrum is assumed to convert from the 
observed count rate or magnitude to the intrinsic flux 
or magnitude. We fit the intrinsic multiband spectrum 
after reddening correction with a power-law function, 
F v oc v a , where v is the pivot frequency of the filter in 
Hz and F v is the flux in erg cm~ 2 s _1 Hz -1 . The same 
a is, in return, used for the assumed spectrum. A cou- 
ple times of iterations are needed in the calculation until 
the fitted a is consistent with the input a. Because we 
are only interested in the continuum component of the 
optical spectrum, data from narrow and medium bands 
are not used due to the possibility that the flux is dom- 
inated by emission lines. The optical spectrum of ULXs 
may be highly variable, thus only quasi-simultaneous ob- 
servations executed within 2 days are used in the fits ex- 
cept for Holmberg II X-l, M81 ULS1, and NGC 5204 
X-l, whose quasi-simultaneous observations only cover 
two bands, and spectra from different epochs are con- 
sistent with a single power-law. A power-law relation is 
adequate to fit all spectra in the sample, except for a few 
sources where excessive flux in the F814W band is seen. 
We compared the properties of the ULX counterparts 
with nearby stars for the three newly identified sources, 
NGC 4559 X-7, M83 IXO 82 and NGC 2403 X-l, in 
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Figure [TJ The ULX counterparts differ from the nearby 
stars in that their spectral energy distributions (SEDs) 
are well fitted with a single power-law, while those of 
the nearby stars show a roll-over at short wavelengths 
suggesting a blackbody-like spectrum with a moderate 
temperature. Also, the nearby stars are less luminous 
than the ULX counterpart. These trends are common to 
the whole sample. The dereddened magnitude in the ob- 
serving band for the compact optical counterpart of each 
ULX from all observations are listed in the Appendix 
(Table EJ. 

2.3. Spectral Classification 

Table [3] gives the absolute visual magnitudes and color 
indices for each quasi-simultaneous observation with at 
least three filters for each source, except for Holmberg II 
X-l, M81 ULS1, and NGC 5204 X-l, for which only two 
bands are available quasi-simultaneously and all their ob- 
servations are listed. The most likely stellar spectral clas- 
sification was inferred for each observation based on the 
magnitude and colors. For most sources, the classifica- 
tions are not consistent with all the colors within one ob- 
servation or the data from different observations suggest 
different spectral types for a given source. The excep- 
tions are IC 342 X-l and NGC 4559 X-7, which appear 
to be consistent with a fixed spectral type. 

As a further check, we compared stellar models for the 
candidate classifications to all of the available data for 
each source and calculated the % 2 deviation of the (non- 
simultaneous) data from the stellar model normalized to 
the average V-magnitude. For NGC 1313 X-2 and M101 
ULX-1, where there are a large number of observations 
available, we chose a subset of the data. The \ 2 an d de- 
grees of freedom (dof) are reported in the table. In most 
cases, the fit to the stellar model is strongly excluded by 
the data. The exceptions are IC 342 X-l and NGC 6946 
X-l. IC 342 X-l appears consistent with an F5I star both 
from the stellar spectral fit and the colors noted above. 
NGC 6946 X-l has a low x 2 because of the large uncer- 
tainty in it reddening correction. The x 2 °f NGC 4559 
X-7, which has a consistent spectral classification based 
on colors, is dominated by variability. 

2.4. Distribution of Spectral Parameters 

Diagrams of My versus a and My versus (B — V)o 
are shown in Figure [2] and the distribution of the opti- 
cal colors is shown in Fig. [3] A value for a is computed 
only when the spectrum has at least three bands, except 
for Holmberg II X-l, M81 ULS1, and NGC 5204 X-l 
whose quasi-simultaneous spectra have only two points. 
The absolute magnitude My is found in a wide range, 
from —7 to —3, for all sources. The spectral index a also 
appears in a wide range from —0.6 to around 2.0, but 
seems to have a distribution peaked between 1 and 2. 
The a of NGC 6946 ULX-1, NGC 5204 X-l and Holm- 
berg IX X-l is consistent with 2.0, which is the slope of 
the Rayleigh- Jeans tail of a blackbody spectrum. The a 
of most sources is less than 2.0, consistent with the fact 
that no thermal emission can produce a spectrum with 
a > 2. It is worth noting that the derived a is sensitive to 
the extinction that is assumed; higher extinction leads to 
larger a. Here, many sources are assumed to have Galac- 
tic extinction, which may be lower than the total extinc- 
tion by 30% or so (see the E g /E n ratios above). Taking 




Fig. 2. — The absolute visual magnitude My versus the spec- 
tral power-law index a (top) or versus (B — V)o (bottom). The 
dotted lines at a = 1/3 and a = 2 are respectively the values ex- 
pected for the intrinsic emission from a steady thin disk and for the 
Rayleigh-Jeans tail of a blackbody. Different points with the same 
color are from different observations of the same source. Sources 
for which the reddening was derived from an associated nebula 
are shown as triangles, while circle denote sources for which only 
Galactic reddening was assumed. The sources are: 1 Holmberg II 
X-l; 2 Holmberg IX X-l; 3 IC 342 X-l; 4 M81 ULS1; 5 M81 X-6; 
6 M83 IXO 82; Y M101 ULX-1; 8 NGC 1313 X-2; 9 NGC 2403 X-l; 
10 NGC 4559 X-7; n NGC 5204 X-l; 12 NGC 5408 X-l; 13 NGC 6946 
ULX-1. 

into account of this factor, it is likely that the peak of a 
distribution is slightly larger and close to 2. The a ex- 
pected for the intrinsic emission from a steady thin disk 
is 1/3, which, however, is inconsistent with most sources 
except NGC 2403 X-l. 

2.5. Long-term Variation 

F555W and F606W are the two most often used filters 
in our sample. Magnitudes in Johnson V were calculated 
from these two bands to examine variability (see Tablc[|)]). 
The observed range in the V-band, U max — U m i n , was 
calculated for each source and is shown in Table |4j 

For M101 ULX-1 and NGC 1313 X-2, multiple quasi- 
simultaneous observations in two bands are available, al- 
lowing us to investigate long-term variability. Not only 
the magnitude, but also the color index varies signifi- 
cantly, see lightcurves of (F555W - F814W) for M101 
ULX-1 and of (F450W - F555W) for NGC 1313 X-2 in 
Figured For M101, we plotted ACS data from 2006 Dec 



HST Photometry of ULX Counterparts 



5 



_ 1 
a> 

§ 

i 10 

8 

6 

4 

2 




(a) ULXs 



(b) LMXBs 



-0.4 -0.2 0.0 0.2 

{B~V) 



0.4 0.6 




.Q 




Fig. 3. — The two upper panels are histograms of the average 
color indices B — V and U — B. The two lower panels are distribu- 
tions of B — V and U — B for optic al counterparts of LMXBs from 
Ivan Pa radiis & McClintock ( 199 5J). The sources indicated in whit e 
include some strange sources in|van Paradiis & McClintock l|1995I V 



23 to 2007 Jan 21. For NGC 1313 X-2, all WFPC2 data 
are plotted. We have checked lightcurves obtained using 
different sky background subtractions and they show the 
same variation pattern, suggesting that the variability 
does not arise from the background. For each ULX, a 
nearby comparison object was chosen to examine instru- 
mental effects. For M101 ULX-1, fitting with a constant 
results in % 2 per degree of freedom (D.o.F.) of 2.6 and 
0.6, respectively for the ULX and the comparison object. 
For NGC 1313 X-2, the color seems to be constant during 
the first three quarters of the total observing window, but 
exhibits large variation in the last 6 observations. Fitting 
with a constant for the last 6 colors, the x 2 /d.o.f. is 3.6 
and 0.8, respectively for the ULX and the comparison 
object. 

Change of the intrinsic emission or the extinction along 
the line of sight could both give rise to the observed 
color variation. Assuming constant intrinsic emission, 
the change of extinction can be directly derived from 
the change of observed magnitude in each band. How- 
ever, for each source, the inferred extinctions at the two 
bands do not agree with each other. Therefore, pure ex- 
tinction variation can be ruled out as the nature of the 
observed variability, and variation of the intrinsic op- 
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Fig. 4.— Color variation of M101 ULX-1 and NGC 1313 X- 
2 (plus) and comparison objects (diamond). For the comparison 
object of M101 ULX-1, its flux on MJD 54096 was affected by 
cosmic rays and thus not plotted. 

tical emission is required. The irregular nature of the 
variability indicates that it cannot be completely caused 
by a periodic phenomena, such as ellipsoidal variations. 
Changes in X-ray illumination, perhaps due to changes 
in the mass accretion rate, could cause the observed vari- 
ab ility. Indeed, for NGC 13 13 X-2, as previously shown 
by llmp iombat oet al.l ((2010), the variability is higher in 
the B-band (F450W) than in the V-band (F555W), sug- 
gesting that the variability arises from X-ray irradiation. 
Unfortunately, there were no suitable, contemporaneous 
X-ray observations that would allow us to study the cor- 
relation between X-ray and optical variability. 

2.6. X-ray to Optical Flux Ratio 

The X-ray to optica l flux ratio, defined following 
iMaccacaro et al.l (|1982D as log(/ x //v) = log/x + 
my/2.5 + 5.37, where fx is the 0.3-3.5 keV observed 
flux in ergs cm -2 s _1 and my is the visual magnitude, 
is calculated for each source. This value can be used 
to distinguish active galactic nuclei (AGN), BL Lac ob- 
jects, clusters of gal axies, normal galax ies, normal stars, 
and X-ray binaries (jStocke et al.l Il99ll ). Another form 
of the X-ray to optical flux ratio as £ = B + 2.5 log Fx 
has been proposed to distinguish between low mass X- 
ray binaries (LMXBs) and high mass X-ray binaries 
(HMXBs), where B is the dereddened B magnitude 
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Fig. 5. — Distributions of the X-ray to optical flux ratio, defined 
as log(/ x //y) or £, for ULXs in the sample, LMXBs, and HMXBs. 
M81 ULS1 is not detected in 2-10 keV in the low state, and thus 
the left arrow indicates its upper limit. 



and £ using the high and low state X-ray fluxes, respec- 
tively. As M81 ULS1 and M101 ULX-1 are very soft and 
radiate most of their fluxes below 2 keV, the 2-10 keV 
flux calculated from models in the literature may be in- 
accurate as model parameters were mainly obtained by 
fitting the spectrum in the low energy band. Therefore, 
we measured Chandra count rates in the 2-7 keV band 
and converted these to fluxes in the 2-10 keV band. This 
conversion was done using the power-law index (but not 
normalization) of the best-fit models from the literature 
and response matrices calculated at the source position. 
M81 ULS1 is not detected in the 2-10 keV band in the 
low state, we thus adopted an upper limit on the flux 
calculated from the measured upper limit on the count 
rate. 

For other sources, their X-ray flux variability is usually 
less than a factor of 10 or even 5, which does not affect the 
X-ray to optical flux ratio significantly. We thus adopted 
a typical model for each source from the liter atur<o and 
calculated the flux given the above energy bands. The 
X-ray fluxes and the X-ray to optical flux ratios are listed 
in Table H 

The distribution of log(/x//v) and £ of ULXs in the 
sample are plotted in Figure [5] The £ distributions of 
LMXBs and HMXBs from Ivan Paradiis fc McCli ntockl 
(1995) are displayed for comparison. We note that the 
3 of the 13 LMXBs and 3 of the 7 HMXBs in the com- 
parison sample are black hole candidates. For HMXBs, 
different types are separately indicated as is in the ref- 
erence. <; of most ULXs in the sample is consistent with 
that of LMXBs, and is statistically higher than the dis- 
tribution of HMXBs. The X-ray to optical ratios of M101 
ULX-1 and M81 ULS1 at the low state are significantly 
lower than those of other ULXs. At the low state, their 
log(/x//y) is around zero, consistent with that of AGN, 
and their £ is lower than that of LMXBs but similar to 
HMXBs. 



and Fx is the 2-10 keV observ ed X-ray flux in /iJy 
van Paradiis fc McClintockl[T995h . 

Due to the highly variable nature of these sources, 
simultaneous X-ray and optical observations may be 
needed to calculate the flux ratio. Holmberg IX X-l 
and NGC 1313 X-2 are the only two sources in our sam- 
ple that simultaneous X-ray and optical observations are 
available. For the Chandra observation of Holmberg IX 
X-l on 2004 February 7, the pileup effect is manageable; 
its X-ray flux is derived by fitting the spectrum with a 
power-law model subject to absorption and pileup. For 
the Chandra observation of NGC 1313 X-2 on 2004 Feb 
22, the source is placed on the chip with a large offset, 
such that the pileup is negligible; an absorbed power-law 
model is used to fit the spectrum. 

M101 ULX-1 and M81 ULS1 are the two most vari- 
able X-ray sources in our sample. The typical X-ray 
luminosity of M101 ULX-1 is about 2 x 10 37 e rg s" 1 
in the low state in 0.3-7 keV (jKong et all 12004ft . or a 
bolom etric luminosity o f 3 x 10 39 erg s _1 in the high 
state (IMukai et al.ll2005l) . For M81 ULS1, the 0.3-8 keV 
flux is 9.3 x 10~ 15 erg s _1 cm -2 in the lo w state an d 
5.3 x 10~ 13 erg s" 1 cnT 2 in the high state (|Liull2008[ ). 
The optical flux does not appear to vary as much as in 
the X-ray band. We thus calculate two sets of log(/x//y) 



3. DISCUSSION 

We have analyzed the largest sample to date of HST 
observations of ULXs with a unique compact optical 
counterpart including 3 new identifications and 10 coun- 
terparts identified in the literature. For each source, we 
measured the broad band optical spectral energy distri- 
bution (SED), its variation, and the X-ray to optical flux 
ratio. In the following, we discuss constraints on the na- 
ture of ULXs and their optical emission derived from our 
results. 

3.1. Origin of the optical emission 

X-ray variability rules out that ULXs are clusters 
of galaxies, normal galaxies, or normal stars. The 
log(/x//y) distribution of ULX in our sample, peaked 
around 3, suggests that most of them are X-ray binaries 
instead of AGN or BL Lac objects, whose log(/x//y) 
ranges from —1 to 1.2 or 0.3 to 1.7, respectively 
(jStocke et al.llT99l . 

3 Holmberg II X-l, IC 342 X-l, and NG C 5204 X-l 
I Feng fc Kaa ret 2009); M81 X-6 (Mizuno el al. 2 0M); M83 I XO 82 
I Stobbart et al.l2006l); NGC 131 3 X-2 JFene fc Ka aret 2007 ); NGC 
2403 X-l IIFene : fc Kaarel l2005h ; NGC 4559 X-7 l lStobbart et al.l 
120061): NGC 5408 X -l [|Kaaret fc Fend I2009T ) : NGC 6946 ULX-1 
UKaiava et al.ll2009T ). 
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Optical light from an X-ray binary can arise from the 
companion star, direct emission from the accretion disk, 
reprocessing of X-rays from the inner disk in the outer 
disk, or some combination of the above. If the optical 
light is dominated by the companion star, then the op- 
tical colors in all observations should be consistent with 
a single spectral type and the optical variability should 
be low. The only variable stars with colors and magni- 
tudes similar to those of most ULX counterparts are /3- 
Cepheid and PV Telesc opii stars; such star have variabil- 
ities of 0.1 mag or less jStankov fc Handl er 2005; Jcffcrv 
2008); even heated companion stars, which would vary in 
an ellipsoidal pattern, are ruled out due to the irregular 
variation of most ULXs. If the optical light is domi- 
nated by direct emission from a multi-color disk black- 
body (MCD), then the optical SED should have a v 1 ^ 
form consistent with extrapolation of the X-ray spectrum 
and there should be significant optical variability. If the 
optical light is dominated by reprocessing of X-ray in the 
outer accretion disk, then the optical SED should have 
a roughl y power-law form with an ind ex between about 
-1 to 2 (|Gierlihski. Done fc Pagd2009ft and, again, there 
should be significant optical variability, correlated with 
the X-ray variability. Also, the X-ray to optical flux ra- 
tio should be similar to those observed from LMXBs, in 
which the optical light is usually dominated by reprocess- 
ing of X-ray in the outer accretion disk when in active 
states. If the companion and disk produce comparable 
optical fluxes, then the variability would be reduced rel- 
ative to that expected from a disk alone and the SED 
would be a combination of a roughly power-law and a 
blackbody form. The observed £, variability, stellar spec- 
trum fitting, and spectral index are used as evaluation 
criteria for different models in Table [5l These properties 
favor the irradiated disk model for most sources. 

The distributions of £, {B-V) , and (U-B) of ULXs 
in the sample appear to be similar to LMXBs, see Fig- 
ure [3] and [5j This does not mean that ULXs actually 
contain low mass companions, but suggests that the op- 
tical spectra of ULXs are dominated by emission from 
the disk instead of the companion star. The fact that 
the whole optical SED of most ULXs cannot be matched 
with a single spectral type, see Table El is further evi- 
dence that the optical emission is not dominated by the 
companion stars. 

Some sources show significant optical variability: M101 
ULX-1, M81 ULS1, NGC 1313 X-2, and NGC 5204 X- 
1. Notably, both sources for which there are more than 
10 observations, M101 ULX-1 and NGC 1313 X-2, show 
significant variability. The non-detection of optical vari- 
ability in the other sources may be due to the limited 
number of observations available. Further, these two 
sources also show color variations that cannot be ex- 
plained as solely due to changes in extinction, there must 
be intrinsic changes in the optical spectra. The detected 
optical variability is 0.2 mag or larger. The only vari- 
able stars with colors and magnitudes similar to those of 
the ULX counterparts are /3-Cepheid and PV Telescopii 
stars . However, such star have variabilities o f 0.1 mag or 
less (jStankov fc Handlerl [20051: IJeffervl[2008h . Thus, the 
variability is inconsistent with a scenario in which the 
companion star dominates the optical light. The vari- 
ability is consistent with the assertion that the optical 



light from ULXs is dominated by emission from the disk. 

The optical spectral index a, see Figure [2] for most 
of the ULXs in our sample lies in the range a = — 1 to 
2. Irra diated disks can produce optica l spectra with such 
slopes ( Gierli hski. Done fc P age 2009). A spectrum with 
a = 2 corresponds to the Rayleigh- Jeans tail of thermal 
emission and would be produced if all of the optically 
emitting regions have temperatures above ~ 2 x 10 4 K 
corresponding to the shortest wavelength of the observed 
optical photons. Spectral indexes between —1 and 2 sug- 
gest that there are optically emitting regions of the disk 
with temperatures down to ~ 6000 K. Variations in the 
optical brightness and spectral slope of the reprocessed 
emission could be produced by variations in the X-ray 
flux from of the central object, the disk radius, or the 
fraction of flux intercepted by the outer disk. 

A B0 lb supergiant companion has been sug gested for 
NGC 5204 X-l based on UV spectroscopy (|Liu et al.1 
2004). We re-examined the data but found that the 
Si in 1299 A line that was used for the classification lies 
close to the strong geocoronal O I 1304A line, and is thus 
suspicious. A background spectrum also shows a similar 
absorption feature. Plus, characteristic absorption lines 
like C iv and Si iv seen in typical B0 lb spectra do not 
appear in the spectrum of NGC 5204 X-l. In the current 
study, the source displays remarkable variation in color 
or spectral index, see source #11 in Figure [21 that could 
not be produced if the optical emission is dominated by 
light from the companion star. Also, the X-ray to op- 
tical flux ratio is consistent with that for a LMXB. We 
therefore argue that the optical/UV emission of NGC 
5204 X-l is dominated by the accretion disk and not by 
a supergiant companion star. 

3.2. Size of the optically emitting region 

If the optical light is due to thermal emission, then 
the size of the emitting region can be estimated from 
the flux density measured at a particular frequency us- 
ing Planck's law if the temperature and distance to the 
source are known. Assuming a simple blackbody model 
for the optical emission, the measurements for Holmberg 
IX X-l suggest an emitting region size of ~ 3 x 10 12 cm 
for a distance of 3.6 Mpc and a blackbody temperature 
of T = 2 x 10 4 K. The size estimate decreases by a fac- 
tor of ~ 5 if the assumed temperature is increased by a 
factor of 10 to T = 2 x 10 5 K. The size depends more 
sensitively on temperature for lower temperatures; it in- 
creases by a factor of 2 if the temperature is decreased 
to T = 1 x 10 4 K and a factor of 9 for T = 5000 K. 
Temperatures lower than 5000 K appear to be excluded 
by the optical spectral shape for most of the sources. 

It is possible to check if the size of the emitting re- 
gion is consistent with a picture in which the opti- 
cal emission is due to reprocessing. Assuming a frac- 
tion rj of the X-ray luminosity Lx is reprocessed, then 
77 = aT A A/Lx where T and A are the temperature and 
area, respectively, of the reprocessing region and a is the 
Stcfan-Boltzmann constant. Inserting T = 2 x 10 4 K, 
A = 1 x 10 24 cm 2 , and Lx = 1 X 10 40 erg s _1 , we find 
77 ~ 0.001. This is similar to values found for the Galac- 
tic black hole X-ray binary XTE J1817-330 when in 
soft sp ectral states similar to those usu ally observed from 
ULXs (jGierlihski. Done fc Pagdl2009|) . Thus, again, the 
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properties of the optical emission appear to be consistent 
with reprocessing. 

The Milky Way black hole binary with the longest 
orbital period, GRS 1915+105, has a separation of 
8 x 10 12 cm and an outer disk radius of 5 x 10 12 cm, 
comparable to the values discussed her e. The 60 day 
X-ray periodicity seen from M82 X-l (jKaaret fc Fend 
120071) is longer than that of GRS 1915+105, suggest- 
ing a larger orbital separation. The binary separation 
of M82 X-l is estim ated to be a few times 10 13 cm 
(jFeng fc Kaaretll2"01Q[ ). In general, large accretion disks, 
of orde r 10 12 cm, would require long orbital periods. Fol- 
lowing [Frank. King, fc Rainel (|2002f ). the disk outer ra- 
dius should be smaller than the circularization radius, 
R c , for matter flowing through the Lagrangian point to- 
wards the black hole. The orbital period is related to R c 



as 



day 



) = ( 



R, 



2.9 x 10 u cm ; 



\3/2 f -3/2 



M 



-1/2 



where Mx is the mass of the accretor and / = (1 + 
g) 4 / 3 [0.5 — 0.227 log q] 4 is a function of the mass ratio 
q = Mq/Mx, and Mq is the companion mass. Requiring 
R > 10 12 cm, for M x = 10M© and q = 0.1, we find 
that P > 11 days. This would be a IM@ companion. 
The period is longer if the companion is heavier, e.g. 
P >26 days for q=0.5. Shorter periods would be possible 
for more massive black holes. 

3.3. Comparison with binary synthesis models 

The My vs. (B — V)o diagram can be compared di- 
rectly with that obtained from binary population syn- 
thesis, in which they assume the optical spectrum com- 
prised of emissi on from both a steady th in disk and 
the companion (Madhusu dhan et al.l 12008ft . Based on 
their models, the most likely color of ULXs is about 
(B — V)o = —0.3 ~ —0.2 independent of the mass 
of the accretor. However, ULXs with an intermediate 
mass black hole are brighter than those with a super- 
Eddington stellar mass black hole of the same luminos- 
ity, with the most probable My w —6 for the former 
and My » —4 for the latter. Our results are consis- 
tent with their model that most ULXs have a (B — V)o 
color about —0.3 to —0.2. The concentration of most 
sources at My rs — 6 would suggest that most ULXs in 
our sample contain intermediate mass black holes. How- 
ever, we note that this result is model dependent. Other 
binary evolution models with different assumptions, i.e. 
that some mass is los t from the system rather t han all 
mass being accreted ()Patruno fc Zampieril 12010ft . could 
lead to different conclusions. 

3.4. NGC 2403 X-l and M83 IXO 82: Intrinsic disk 
emission? 

NGC 2403 X-l has an optical spectral index con- 
sistent with that expected, a — 1/3, for the intrin- 
sic emission of a standard multicolor disk. Assuming 
a = 1/3, its optical spectrum can be written as F v = 
6 x 10~ 35 iW 3 erg cm~ 2 s _1 Hz -1 . A face-on standard 
disk at 3.2 Mpc with an inner temperature fcTj n = 1 keV 
and a bolometric luminosity Lboi = 3 x 10 40 erg s _1 , or 
with kT m = 0.1 keV and Lboi = 10 39 erg s _1 , can pro- 
duce an optical tail consistent with the observed spec- 
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Fig. 6. — Variable optical emission from M101 ULX-1. The dot- 
ted line indicates a power-law model fitted to the short wavelength 
points with a = 1.52 ± 0.12 for the data from 1994 Apr 08. 



trum. The measured X-ray lumi nosity for the source 
ranges from 3 to 9 x 10 39 erg s" 1 (jFeng : fc Kaareti[2009l) 
which is in the required range. Thus, the optical emission 
for this source may be dominated by intrinsic emission 
from an accretion disk. However, simultaneous X-ray 
and optical/UV observations are needed to test this hy- 
pothesis. 

The optical spectrum of M83 IXO 82 could be consis- 
tent with the a = 1/3 expected for intrinsic emission of 
a standard multicolor disk if our assumed reddening is 
incorrect. We calculated the reddening assuming only 
Galactic absorption of E(B — V) = 0.066. Additional 
intrinsic absorption of E(B — V) = 0.111 would make 
a consistent with 1/3. This source shows an absolute 
visual magnitude, My ~ —3, and an X-ray luminosity, 
Lj~ 2 x 10 39 erg s -1 ()Feng fc Kaaretl [20051 ). similar to 
that of NGC 2403 X-l. Thus, this source is also a good 
candidate to have optical emission dominated by intrin- 
sic emission from an accretion disk. 

3.5. M101 ULX-1 and M81 VLSI 

M101 ULX-1 and M81 ULS1 are exceptional cases ow- 
ing to their strong X-ray variability and deviation from 
the main population in the X-ray to optical flux ratio dis- 
tributions. Especially when in the low state, their X-ray 
to optical flux ratios are consistent with that of AGN. 
However, detection of variability on time scales o f 10 3 s 
from both sources (jMukai et all 120051: iLiul 12008ft . more 
rapid than the X-ray variability seen from AGN, makes 
an AGN identification unlikely. The X-ray to optical flux 
ratio of the two sources is consistent with that of HMXBs. 
Many Galactic HXMBs, especially with a neutron star 
accretor, are transient sources with huge variability due 
to uneven wind- fed accretion in an elliptical orbit. The 
high luminosity and moderate absorption column density 
of these ULXs suggest that they are Roche-lobe overflow 
system, where circular orbits are often presumed. How- 
ever, the large long-term variability of these two systems 
may suggest eccentric orbits. 

The optical spectral properties of the two sources are 
different. The SED of M101 ULX-1 at 3 epochs is shown 
in Fig. [6] There is variability at every waveband where 
there are multiple observations and the variations are 
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Fig. 7. — Optical spectra of IC 342 X-l from two sets of quasi- 
simultaneous observations. The lines are the power-law and black- 
body model fitted to the Dec 2005 data. The arrow indicates the 3<r 
upper limit of the flux at F330W where the source is not detected. 



stronger at long wavelengths. The variability could be 
explained in a disk irradiation as due to changes in the 
outer disk radius. If the blue wing of the spectrum is 
due to companion ligh t, then a ve ry hot star would be 
required. We note that ILiu I ()2009D fitted the same data 
using the sum of a Wolf-Rayet star, an irradiated disk 
described by an a = —1 spectrum, and an additional 
hot component ascribed to the X-ray heated side of the 
compan ion. We note that more compl ex irradiated disk 
models (|Gierlihski. Done fc Page 20( M) can produce the 
extra hot component needed bv ILiu I (|2009T ) in the inner 
regions of the disk suffice without need of an additional 
physical surface. 

In contrast, the optical spectrum of M81 ULSI is signif- 
icantly shallower. At least for one observation, the spec- 
trum could be consistent with a = 1/3 as expected for 
intrinsic emission of a standard multicolor disk if there is 
additional reddening beyond t he Galacti c value assumed 
in our analysis. We note that ILiu et ail (|2008D modeled 
the blue part of the optical spectrum as due to direct 
disk emission. However, this modeling was based on si- 
multaneous fitting of HST observations obtained over a 
period of 4 years. Given the highly variable nature of 
the source, it is important to better constrain the spec- 
tral behavior with simultaneous observations extending 
into the UV. 

3.6. IC 342 X-l: A single blackbody? 

The spectrum from IC 342 X-l on December 2005 
can be equally well fitted by a power-law spectrum with 
a = —0.62 ± 0.18, or a blackbody spectrum with a tem- 
perature of 7100±400 K, see Figure[?l However, the spec- 
trum from September 2005 shows a cutoff in the F330W 
band, inconsistent with the power-law. Fitting to the 
combined data favors the blackbody model. 

Such blackbody emission could arise from a companion 
star. There is no evidence of optical variability and a 
spectral classification of F5 lb or lab fits all the data 
available at both observation epochs. If the reddening 
is less than the nebular value we have assumed, then 
the star would be a later spectral type, i.e. F8-G0 lb 
(|Feng fe Kaaretl I2008D . Such blackbody emission could 
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Fig. 8. — The spectrum of Holmberg IX X-l shows a red excesses 
in the F814W band. The dotted line indicates a power-law model 
fitted to the short wavelength points, with a = 2.10 ± 0.23. 

also arise from irradiation in a very large disk. For IC 342 
X-l, the blackbody fitting suggests a surface area of 6 x 
10 26 cm 2 , or a reflector radius of 10 13 cm, which is about 
an order of magnitude larger than the value inferred for 
the majority of ULXs. Observations at multiple epochs 
to test for optical variability could distinguish between 
these two possibilities. 

3.7. Holmberg IX X-l: Red excess 

The optical spectrum of Holmberg IX X-l is well fit- 
ted at short wavelengths by a single power-law with a 
spectral index a = 2.10 ± 0.23. This may represent 
the Rayleigh- Jeans tail of reprocessed disk emission with 
temperatures above about 2 x 10 4 K. At the red end, the 
spectrum shows an excess at the F814W band above the 
extrapolation of the power-law spectrum at short wave- 
lengths (Figure [5J . This indicates some source of emis- 
sion at cooler temperatures and could be due either to 
the companion star or some structure in the disk at large 
radii. New observations with multiple bands at several 
epochs could distinguish between these two possibilities. 



3.8. NGC 5204 X-l: 



not associated with nearby 
nebulae? 



NGC 5204 X-l appears o n the sky near an optical 
nebulae (jPakull et al.l [2002) wi th a measured nebular 
extinction of ' E(B - V) = 0.4 (jAbolmasov et al.l 120071: 
iRussell et alj 12010). The Galactic extinction along the 
line of sight of the source is E(B — V) = 0.013, which 
is adopted for reddening correction in this paper. The 
source spectrum is found to vary and the power-law index 
a is around 2.0, which is the largest value theoretically 
allowed. Adopting the nebular extinction would lead to 
an unphysical steep spectral slope, with a > 2.8. The 
X-ray absorption column density of the source is about 
(0.6 - 1.0) x 10 21 cm" 2 (jFeng fc Kaaretl [2009]) . which 
translates to an extinctio n of E{B — V) = 0.11 — 0.22 
(|Predehl fc Schmitt]|1995[ ) al lowing for the range of dust- 
to-gas rati os discussed in iCardelli. Clavton fc Mathisl 
((19891) and IPredehl fc Schmitti (119951). This is signifi- 
cantly lower than the nebular value. Therefore, we sug- 
gest that the X-ray source may be not associated with 
the nebulae. 
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4. CONCLUSIONS 

• The dominant component of optical emission in most 
ULXs is produced via X-ray irradiation of the outer 
accretion disk, similar to the case of LMXBs. 

• As a corollary, one cannot use broad-band optical col- 
ors to infer the companion spectral type. 

• Large optical emitting regions are required, suggesting 
relatively large orbital separations and long orbital pe- 
riods. 

• The optical spectrum and X-ray luminosity of NGC 
2403 X-l are consistent with both arising from intrinsic 
emission from a standard accretion disk. M83 IXO 
82 is a somewhat weaker candidate for intrinsic disk 
emission extending into the optical. 

• M101 ULX-1 and M81 ULS1 are unique in high X-ray 
variability and low X-ray to optical flux ratios. These 
systems may have eccentric orbits. 

• IC 342 X-l has a spectrum quite different from that of 
most ULXs and the optical emission may be dominated 



by an F5 or somewhat later supergiant companion star. 
• NGC 5204 X-l seems not to be physically related to the 
nearby optical nebulae due to inconsistent extinction. 
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TABLE 3 

Absolute visual magnitudes, colors, and spectral classifications of the ULXs in the sample 



ULX 


date 


M v 


U-B 


B-V 


V — I 


V - R Sp. 


X 2 /dof 


Holmbcrg II X-l 


2006-12-30 


-6.00 






-0.02 


B6-B7 lb-lab 


B0 I, 112.4/8; B0. 5 I, 70.6/8; 




onn7 i r\ no 


—5.93 




— 0.24 




BU-B1 lb 


DP T TO K / Q 

Bo 1, fO.O/O 




nnnT 1 ( \ ( \v 

2007-10-05 


—5.94 




—0.18 




B1-B2 lb 




1 n t\C\ 

2007-10-09 


—5.94 




—0.27 




An nn tu 

(J9-B0 lb 




TT 1 1 TV V 1 

Holmbcrg IX X-l 


2004-02-07 


—5.88 


— 1.33 


—0.42 


—0.04 


nn r n^ni 

05V, 06111 


/""\ o ~\ t ooo o I A f \ v AT" or*o o//i 

03 V, 2AI.A/A; U5 V, 253.3/4 


IC 342 X-l 


2005-09-02 


-5.94 


> -0.77 


0.37 




F5 lb-lab 


F5 I, 10.7/5 




onncr i n i o 

2U05-12-18 


—5.95 




0.30 


0.60 


U.oo r 5 lb-lab 


M81 ULS1 


2006-03-21 


-6.26 




0.11 




A5-A7 lab 


F2 I, 424.1/6; F5 I, 418.3/6; 




nr\r\a no 07 
zUUu-Uo-z / 


—5.99 




0.22 




• • • r 2 lb-lab 


ttq t 7oc o la.. (~*c\ T i ccncc i la. 
ro 1, foo.o/o; Kj\j 1, loUo.l/o 


iv toi "v a 


i nn c m n 
1995-01-01 


—3.91 


— 1.51 


0.06 


—0.05 


n r\r~ i~> cr A o tt 

— 0.0d Bo-Ao 11 


B2 II, 78.6/14; B5 II, 61.6/14 




2001-06-04 


—4.01 




—0.13 


0.07* 


B7-B8 11 


moo tv/~\ oo 

Moo 1X(J o2 


o r\ ri no o tr 

zUUo-U2-z5 


—3.11 


—0.96 


0.26 




AO TTin TT T"> 1 \T 

A2-r 2 11, Bl V 


rin tt i u o / o rio tt i"i cr/o 

1*0 11, lo.z/z; b2 11, ll.o/z; 
Bl V, 21.7/2 


TV/Tim T TT V 1 

M101 ULX-1 


i / \ * i i no nn 

1994-03-22 


—6.39 




—0.11 


0.20 


B5 lab, A3 lab 


i r ~\ x ^ino/n o TTT O/in o /n 

05 1, 412.8/9; U8 111, 342.3/9; 




1994-04-08 


—6.15 


— 1.20 


—0.28 


0.10* 


Ar TTT /~~\ O TT 

05 III, 08 II 


B5 I, 320.1/9 




o r\ c\ n 1 1 1 r 

2002-11-15 


—6.00 




—0.31 


—0.01* 


s~\ f TTT ^\C> j^An TT T~» >~> Tl 

05 111, 08-09 11, B7 lb 


JNGC 1313 X-2 


2003-11-22 


—4.61 


— 1.14 


—0.26 


—0.09* 


(~\ / \T j T) n T> 1 TTT 

U9V, B0-B1 111 


j~\r\ \7 <~?r* a /o T) n TTT rrx 1-7 1 

09 V, 76.4/3; B0 ill, 62.7/3; 
Bl III, 153.1/3 


NGC 2403 X-l 


2005-10-17 


-2.90 


-1.01 


0.07* 




B0-B2 V 


B0 V, 33.6/2; Bl V, 28.5/2; 
B2 V, 33.3/2 


NGC 4559 X-7 


2001-05-25 


-6.98 




-0.17 


-0.11 


B2-B5 la 


B5 I, 26.0/6 




2005-03-08 


-7.10 




-0.12 


-0.06 


B3-B5 la 


NGC 5204 X-l 


2001-05-28 


-5.77 






-0.52 


05V 


09 III, 33.8/9; B0 I, 25.7/9; 




2008-08-08 


-5.57 




-0.35 




05V, 08III 


05 V, 707.5/9; B2 II, 1060.8/9 




2008-08-10 


-5.62 




-0.23 




B2 Il-Ib 




2008-08-13 


-5.69 




-0.20 




B3 Il-Ib 




NGC 5408 X-l 


2009-04-04 


-6.27 


-1.11 


-0.28 


-0.19* 


09 lb 


08 I, 29.8/5; 09.5 I, 37.1/5 


NGC 6946 ULX-1 


1996-01-27 


-7.26 




-0.07 




B6-B7 la 


05 I, 8.3/5; 08 I, 7.3/5; 




2001-06-08 


-7.27 




-0.42* 


-0.21 


B2 la 


B0 I, 4.3/5; B5 I, 6.0/5 



Note. — * Color that does not agree with the best-fit spectral type. 



TABLE 4 

Optical and X-ray Properties of the ULXs in the sample 



ULX 


date 


My 


a 


X 2 /dof 


fx 




Fx 




log(/x//v) 




Vmax Vrnin 


# obs 


Holmbcrg II X-l* 


2006- 12-30 

2007- 10-03 
2007-10-05 
2007-10-09 


-6.00 
-5.93 
-5.94 
-5.94 


1.10 ±0.09 
1.41 ±0.22 
1.17 ±0.21 
1.52 ±0.19 




4.08E- 


12 


1.09E- 


-01 


2.7 


18.9 


0.07 ±0.03 


4 


Holmberg IX X-l 


2004-02-07 


-5.88 


1.79 ±0.19 


44.9/2 


4.92E- 


12 


3.62E- 


-01 


3.2 


20.4 


0.09 ±0.04 


2 


IC 342 X-l 


2005-12-18 


-5.95 


-0.62 ±0.18 


2.0/2 


1.79E- 


12 


1.78E- 


-01 


3.3 


20.1 


0.01 ±0.13 


2 


M81 ULS1* 


2006-03-21 


-6.26 


0.02 ±0.13 




5.43E- 


15 l 


< 3.51E-05' 


-0.1 ! 


< 10.7' 


0.27 ±0.06 


2 




2006-03-27 


-5.99 


-0.39 ±0.13 




2.03E- 


13 h 


5.21E- 


-05 h 


lA h 


11. l' 1 






M81 X-6 


1995-01-31 
2001-06-04 


-3.91 
-4.01 


1.19 ±0.18 
0.94 ±0.14 


3.8/3 
0.1/1 


1.55E- 


12 


1.02E- 


-01 


3.2 


21.3 


0.18 ±0.09 


5 


M83 IXO 82 


2006-02-25 


-3.11 


-0.43 ± 0.35 


0.4/1 


3.18E- 


13 


8.88E- 


-03 


3.0 


20.4 






M101 ULX-1 


1994-03-22 


-6.39 


0.62 ±0.10 


1.0/1 


2.64E- 


15 ; 


1.64E- 


-05 1 


0.2' 


10.9' 


0.67 ±0.03 


27 




1994-04-08 


-6.15 


1.20 ±0.08 


10.8/2 


3.39E- 


13 h 


3.37E- 


-04 h 


2.3 h 


14.2'' 








2002-11-15 


-6.00 


1.41 ±0.07 


7.4/1 


















NGC 1313 X-2 


2003-11-22 


-4.61 


1.37 ±0.14 


3.5/2 


1.72E- 


12 


1.17E- 


-01 


3.1 


20.5 


0.20 ±0.04 


22 


NGC 2403 X-l 


2005-10-17 


-2.90 


0.34 ±0.20 


0.8/1 


8.33E- 


13 


3.81E- 


-02 


3.2 


21.1 






NGC 4559 X-7 


2001-05-25 


-6.98 


1.25 ±0.07 


1.2/1 


5.89E- 


13 


2.24E- 


-02 


2.3 


18.7 


0.11 ±0.04 


3 




2005-03-08 


-7.10 


1.06 ±0.10 


0.5/2 


















NGC 5204 X-l* 


2001- 05-28 

2002- 10-28 
2008-08-08 
2008-08-10 
2008-08-13 


-5.77 

-5.57 
-5.62 
-5.69 


2.42 ± 0.39 
0.92 ±0.11 
1.87 ±0.27 
1.38 ±0.46 
1.26 ±0.28 


1.58E- 


12 


4.95E- 


-02 


2.6 


19.2 


0.12 ±0.03 


4 


NGC 5408 X-l 


2009-04-04 


-6.27 


1.41 ±0.14 


3.6/2 


3.32E- 


12 


2.33E- 


-02 


2.8 


17.8 


0.03 ±0.05 


2 


NGC 6946 ULX-1 


2001-06-08 


-7.27 


1 SS+ ' 48 


0.5/1 


4.77E- 


13 


2.09E- 


-02 


2.2 


16.8 


0.01 ±0.17 


2 


Note. — * Sources that have quasi 


simultaneous observations available 


only at two bands; X 


ray low state; ^X-ray 


high state; 


a is the power- 


law index 



of the optical spectrum; x 2 /dof is derived from fitting with a power-law model; /x is the observed 0.3-3.5 keV flux in erg cm -2 s _1 ; i*x 1S the 2-10 keV 
observed X-ray flux in fiJy: V ma x — Vmin indicates the largest variation in the V band; # obs is the number of observations for which a V magnitude was 
measured. 
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TABLE 5 
Evaluation of different models 



ULX 






Properties 






Models 






Variability 


Stellar spectrum 


Spectral index 


Star 


Irradiated disk 


Direct disk 


Holmberg II X-l 


L 


M 


B 


K 




+ 




noirriDcrg ia yv-± 


T 

J_j 


K I 


r 
i_> 


r» 




i 

-r 




IC 342 X-l 


L 


N 


G 


R 


++ 


± 




M81 ULS1 


H 


V 


B 


R 




+ 




M81 X-6 


L 


M 


B 


R 




+ 




M83 IXO 82 


L 




B 


R 




+ 




M101 ULX-1 


H 


V 


B 


R 




+ 




NGC 1313 X-2 


L 


V 


B 


R 




+ + 




NGC 2403 X-l 


L 




B 


D 







+ 


NGC 4559 X-7 


L 


M 


B 


R 




+ 




NGC 5204 X-l 


L 


V 


B 


R 




++ 




NGC 5408 X-l 


H 


N 


B 


R 









NGC 6946 ULX-1 


H 


N 


G 


R 


++ 








Note. — £: these sources whose £ are larger than 18 are marked as L. otherwise, are marked as H; Variability: V, M and N are used for these 
sources whose V max — V m i n ^ 4 error, 4 error > V max — V m i n ^ 2 error and Vmax — V m i n < error, respectively; Stellar spectrum: using x 2 /dof 
in Table [3] to calculate the probability P that a random variable x 2 ' 1S l ess than or equal to observed \ 2 ■ The stellar spectrum fitting is bad (B) 
if P > 99%, otherwise, is good (G); Spectral index: D is used if a is consistent with 1/3 within errors, otherwise, R is used. Star: H — h if stellar 
spectrum fitting is good and variability is N, — if stellar spectrum fitting is bad or variability is V; Irradiated disk: H — h if £ — L Sz spectral index 
— R Sz variability — V, + if £ — L Sz spectral index — R Sz variability / V or ^ - H & spectral index — R & variability — V, if £ — L Sz spectral 
index ^ R Sz variability ^ V or £ — H 8z spectral index — R Sz variability ^ V; Direct disk: + if spectral index — D, — if spectral index — R. The 
rate of how strongly the model is recommended should be based on the decreasing from ++ to — . 

TABLE 6 

HST OBSERVATIONS AND MAGNITUDES OF ULXS IN OUR SAMPLE 



ULX 


date 


datasct 


instrument and filter 


exposure (s) 


m Filter 


my 



Holmberg II X-l 



2006-01-28 
2006-12-30 

2006- 12-30 

2007- 10-03 
2007-10-03 
2007-10-05 
2007-10-05 
2007-10-09 
2007-10-09 



j9dr03010 
j9cmz7010 
j9cmz7020 
u9zs240 
u9zs240 
u9zs250 
u9zs250 
u9zs260 
u9zs260 



1 


2 


m 


3 


1 


m 


1 


2 


m 


3 


1 


m 


1 


2 


m 


3 


1 


in 



ACS/WFC /F814W 600 21.450 ±0.029 

ACS/WFC /F555W 4660 21.403 ±0.030 21.418 ± 0.030 

ACS/WFC /F814W 4660 21.451 ±0.019 

WFPC2/PC1/F450W 500/500 21.287 ±0.018 

WFPC2/PC1/F555W 500/500 21.471 ±0.012 21.489 ±0.012 

WFPC2/PC1/F450W 500/500 21.328 ±0.015 

WFPC2/PC1/F555W 500/500 21.465 ±0.012 21.478 ±0.012 

WFPC2/PC1/F450W 500/500 21.259 ±0.014 

WFPC2/PC1/F555W 500/500 21.464 ±0.012 21.484 ± 0.012 



Holmberg IX X-l 


2004-02-07 


j8ol03010 


ACS/WFC /F814W 


1160 


21.885 ±0.042 






2004-02-07 


j8ol03040 


ACS/WFC /F435W 


2520 


21.461 ±0.021 






2004-02-07 


j8ola3010 


ACS/HRC/F330W 


2760 


19.616 ±0.024 






2004-02-07 


j8ol03030 


ACS/WFC/F555W 


1160 


21.868 ±0.028 


21.901 ±0.028 




2004-03-25 


j8ol07010 


ACS/WFC /F555W 


2400 


21.960 ±0.021 


21.992 ± 0.021 


IC 342 X-l 


2005-09-02 


j9h8b3010 


ACS/HRC/F330W 


2900 


> 20.868 






2005-09-02 


j9h813010 


ACS/WFC /F435W 


1248 


22.019 ±0.132 






2005-09-02 


j9h813020 


ACS/WFC /F606W 


1248 


21.521 ±0.079 


21.651 ±0.100 




2005-12-16 


j9jn[0|a|b|c]1010 


ACS/WFC /F625W 


920/920/920/920 


21.293 ±0.071 






2005-12-18 


j9jn02020 


ACS/WFC /F814W 


1080 


21.022 ±0.094 






2005-12-18 


j9jn02040 


ACS/WFC /F435W 


1800 


21.942 ±0.094 






2005-12-18 


j9jn02010 


ACS/WFC /F555W 


1080 


21.667 ±0.089 


21.640 ±0.089 



M81 ULS1 2003-09-18 j8mxl8akq ACS/WFC/F814W 120 20.965 ± 0.050 

2004- 09-15 j901a7010 ACS/WFC/F814W 1650 21.019 ±0.051 

2005- 12-09 u9el200[3|4]m WFPC2/WF4/F300W 500/500 20.369 ± 0.075 

2006- 03-21 j9ell5020 ACS/WFC/F606W 1200 21.460 ±0.034 
2006-03-21 j9ell5010 ACS/WFC/F435W 1200 21.648 ±0.020 
2006-03-27 j9el28a9q ACS/WFC /F435W 465 22.026 ± 0.024 
2006-03-27 j9el28aaq ACS/WFC /F606W 470 21.699 ±0.032 



21.541 ±0.044 



21.811 ±0.041 



M81 X-6 



1995 


01 


31 


1995 


01 


31 


1995 


01 


31 


1995 


01 


31 


1995 


01 


31 


2001 


06 


04 


2001 


06 


04 


2001 


06 


04 


2004 


09 


15 


2004 


09 


25 


2004 


09 


25 


2006 


03 


22 


2006 


03 


22 



u2mh010[b|c]t 

u2mh010[d|e]t 

u2mh010[5|6 7]t 

u2mh010[f|g]t 

u2mh010[8|9 a]t 

u6eh010[5r|6m|7r|8r] 

u6eh010[l|2|3|4]r 

u6eh010t9|a|b|c]r 

j90108010 

u8ZO010[l|3]m 

u8ZO010[2|4]m 

j9ell8010 

j9cll8020 



WFPC2/WF4/F555W 

WFPC2/WF4/F675W 

WFPC2/WF4/F336W 

WFPC2 /WF4 /F814W 

WFPC2 /WF4 /F439W 

WFPC2 /PC1/F450W 

WFPC2/PC1/F555W 

WFPC2/PC1/F814W 

ACS/WFC/F814W 

WFPC2/WF3/F675W 

WFPC2/WF3/F555W 

ACS/WFC /F435W 

ACS/WFC/F606W 



300/600 
300/600 
400/160/600 

300/600 
400/400/400 
500/500/500/500 
500/500/500/500 
500/500/500/500 
1650 
600/600 
500/600 
1200 
1200 



23.880 
23.958 
21.939 
23.955 
23.982 
23.674 
23.782 
23.737 
23.560 
23.618 
23.879 
23.722 
23.664 



± 0.060 
±0.107 
±0.142 
±0.109 
±0.166 
± 0.043 
±0.060 
±0.077 
±0.099 
±0.074 
±0.076 
±0.029 
±0.054 



23.894 ± 0.060 

23.791 ± 0.060 

23.872 ± 0.076 
23.717 ±0.068 



HST Photometry of ULX Counterparts 
TABLE 6 — Continued 
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ULX date 


datasct 


instrument and filter 


exposure (s) 


m-Filtcr 


my 


onnfi o^ 97 


j9cla8acq 






ACS /WFC /FfinfiW 




93 7sq -4- n 


93 &9Q -4- n nfis 




j9cla8acq 






ACS /WFC /F435W 


400 


9Q 7S1 -U n 043 




M83 1X082 2006-02-25 


j9h811020 




ACS/WFC /F606W 


1000 


25.132 ±0.107 




2006-02-25 


j9h811010 




ACS/WFC/F435W 


1000 


25.512 ±0.080 




2006-02-25 


j9h8al010 




ACS/HRC /F330W 


2568 


24.156 ±0.263 




M101 ULX1 1994-03-22 


u2780101t 




WFPC2 /WF4 /F555W 


1200 


22.895 ±0.026 


22.898 ± 0.026 


1994-03-22 


u2780102t 




WFPC2 /WF4 /F814W 


1000 


22.712 ±0.041 




1994-03-23 


u2780103t 




WFPC2/WF4/F439W 


1000 


22.813 ±0.070 




1994-03-30 


u2780201t 




WFPC2/WF4/F555W 


1200 


23.022 ± 0.032 


23.040 ± 0.032 


1994-04-08 


u2780302t 




WFPC2/WF4/F814W 


1200 


23.044 ± 0.049 




1994-04-08 


u2780301t 




WFPC2/WF4/F555W 


1200 


23.120 ±0.031 


23.133 ±0.031 


1994-04-08 


u2780303t 




WFPC2/WF4/F439W 


1200 


22.884 ± 0.068 




1994-04-08 


u278030[4|5]t 


WFPC2/WF4/F336W 


1200/1200 


21.143 ±0.058 




1994-04-11 


u278040[lp|2t] 


WFPC2/WF4/F555W 


1200/120 


23.009 ± 0.059 


23.027 ± 0.059 


1994-04-16 


u2780601t 




WFPC2/WF4/F555W 


1200 


23.108 ±0.031 


23.126 ±0.031 


1994-04-20 


u2780701t 




WFPC2/WF4/F555W 


1200 


23.160 ±0.032 


23.178 ±0.032 


1994-04-24 


u2780801t 




WFPC2/WF4/F555W 


1200 


23.192 ±0.040 


23.210 ± 0.040 


1994-04-28 


u2780901t 




WFPC2/WF4/F555W 


1200 


23.167 ±0.028 


23.185 ±0.028 


1994-05-03 


u2780A01t 




WFPC2/WF4/F555W 


1200 


23.222 ± 0.033 


23.240 ± 0.033 


1994-05-05 


u2780501t 




WFPC2/WF4/F555W 


1200 


23.249 ± 0.031 


23.267 ± 0.031 


1994-05-05 


u278050[2|3]t 


WFPC2/WF4/F814W 


1200/120 


23.252 ±0.070 




1994-05-10 


u2780[b|c]01t 


WFPC2/WF4/F555W 


1200/1200 


23.256 ± 0.022 


23.274 ± 0.022 


1994-05-10 


u2780b02t 




WFPC2/WF4/F814W 


1200 


23.316 ±0.053 




1995-03-22 


u2780d01t 




WFPC2/WF4/F555W 


1200 


23.237 ± 0.034 


23.255 ±0.034 


1995-04-17 


u2ms030[l|2]t 


WFPC2/WF4/F555W 


500/500 


23.208 ± 0.039 


23.226 ± 0.039 


2002-11-15 


j8d601031 






ACS /WFC /F814W 


720 


23.315 ±0.047 




2002-11-15 


j8d601021 






ACS/WFC /F555W 


720 


23.261 ±0.030 


23.283 ±0.030 


2002-11-15 


j8d601011 






ACS /WFC /F435W 


900 


22.951 ±0.018 




2006-12-23 


j9o401020 






ACS/WFC /F555W 


1330 


23.177 ±0.014 


23.199 ±0.014 


2006-12-23 


j9o401010 






ACS /WFC /F814W 


724 


23.146 ±0.033 




2006-12-24 


j9o402020 






ACS/WFC /F555W 


1330 


23.167 ±0.015 


23.190 ±0.015 


2006-12-24 


j9o402010 






ACS/WFC /F814W 


724 


23.116 ±0.029 




2006-12-25 


j9o403020 






ACS/WFC /F555W 


1330 


23.205 ± 0.015 


23.227 ±0.015 


2006-12-25 


j9o403010 






ACS /WFC /F814W 


724 


23.144 ±0.027 




2006-12-26 


j9o404020 






ACS /WFC /F555W 


1330 


23.181 ±0.014 


23.204 ±0.014 


2006-12-26 


j9o404010 






ACS/WFC /F814W 


724 


23.183 ±0.032 




2006-12-28 


j9o405020 






ACS/WFC /F555W 


1330 


23.156 ±0.015 


23.178 ±0.015 


2006-12-28 


j9o405010 






ACS /WFC /F814W 


724 


23.227 ±0.030 




2006-12-30 


j9o406020 






ACS/WFC /F555W 


1330 


23.201 ±0.016 


23.224 ±0.016 


2006-12-30 


j9o406010 






ACS /WFC /F814W 


724 


23.119 ±0.030 




2007-01-01 


j9o407020 






ACS/WFC /F555W 


1330 


23.243 ± 0.015 


23.266 ±0.015 


2007-01-01 


j9o407010 






ACS /WFC /F814W 


724 


23.216 ±0.030 




2007-01-04 


j9o408020 






ACS /WFC /F555W 


1330 


23.249 ±0.016 


23.272 ±0.016 


2007-01-04 


j9o408010 






ACS /WFC /F814W 


724 


23.221 ± 0.030 




2007-01-07 


j9o409020 






ACS/WFC /F555W 


1330 


23.311 ±0.017 


23.334 ±0.017 


zUU ( -Ul-U ( 


j9o409010 






Ape /XAfPr^ /T?S1 AWT 
0/ Wr L'/r 514 VV 


70/1 


00 1 v7 -U n riQi 




2007-01-11 


j9o410020 






ACS /WFC /F555W 


1330 


23.179 ±0.016 


23.202 ±0.016 


2007-01-11 


j9o410010 






ACS /WFC /F814W 


724 


23.069 ±0.028 




2007-01-17 


j9o411020 






ACS/WFC /F555W 


1330 


23.009 ± 0.012 


23.031 ±0.012 


2007-01-17 


j9o411010 






ACS /WFC /F814W 


724 


22.928 ± 0.025 




2007-01-21 


j9o412020 






ACS /WFC /F555W 


1330 


23.081 ±0.014 


23.104 ±0.014 


2007-01-21 


j9o412010 






ACS /WFC /F814W 


724 


23.054 ± 0.031 




2008-01-05 


ua30640[3|4]m 


WFPC2/WF3/F555W 


500/500 


22.642 ± 0.024 


22.660 ± 0.024 


2008-01-05 


ua30640[l 2]m 


WFPC2/WF3/F814W 


350 


22.616 ±0.048 




NGC 1313 X-2 2003-11-22 


j8ola2010 






ACS/HRC/F330W 


2760 


21.317 ±0.029 




2003-11-22 


j8ol02040 






ACS/WFC/F435W 


2520 


22.955 ±0.012 




2003-11-22 


j8ol02010 






ACS /WFC /F814W 


1160 


23.339 ± 0.023 




2003-11-22 


j8ol02030 






ACS /WFC /F555W 


1160 


23.211 ±0.017 


23.233 ±0.017 


2004-02-22 


j8ol06010 






ACS/WFC /F555W 


2400 


23.086 ± 0.015 


23.108 ±0.015 


2008-05-21 


u9zs010 


3 


4 


111 


WFPC2/PC1/F555W 


400/700 


23.057 ± 0.028 


23.074 ± 0.028 


2008-05-21 


u9zs010 


1 


2 


m 


WFPC2/PC1/F450W 


500/500 


22.982 ± 0.044 




2008-05-22 


u9zs020 


3 


4 


m 


WFPC2/PC1/F555W 


400/700 


23.149 ±0.033 


23.166 ±0.033 


2008-05-22 


u9zs020 


1 


2 


111 


WFPC2/PC1/F450W 


500/500 


22.881 ±0.054 




2008-05-23 


u9zs030 


1 


2 


111 


WFPC2/PC1/F450W 


500/500 


22.863 ± 0.046 




2008-05-23 


u9zs030 


3 


4 


111 


WFPC2/PC1/F555W 


400/700 


23.115 ±0.034 


23.132 ±0.034 


2008-05-24 


u9zs040 


3 


4 


111 


WFPC2/PC1/F555W 


400/700 


23.153 ±0.030 


23.170 ±0.030 


2008-05-24 


u9zs040 


1 


2 


111 


WFPC2/PC1/F450W 


500/500 


23.029 ± 0.044 




2008-05-25 


u9zs050 


3 


4 


111 


WFPC2/PC1/F555W 


400/700 


23.212 ±0.031 


23.229 ±0.031 


2008-05-25 


u9zs050 


1 


2 


m 


WFPC2/PC1/F450W 


500/500 


23.056 ± 0.053 




2008-05-26 


u9zs060 


1 


2 


111 


WFPC2/PC1/F450W 


500/500 


23.030 ± 0.049 




2008-05-26 


u9zs060 


3 


4 


m 


WFPC2/PC1/F555W 


400/700 


23.254 ± 0.038 


23.271 ±0.038 


2008-05-27 


u9zs070 


1 


2 


111 


WFPC2/PC1/F450W 


500/500 


23.073 ± 0.046 




2008-05-27 


u9zs070 


3 


4 


m 


WFPC2/PC1/F555W 


400/700 


23.188 ±0.031 


23.205 ± 0.031 
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Tao et al. 
TABLE 6 — Continued 



ULX date 


datasct 


instrument and filter 


exposure (s) 


Tft-Filter 


my 


2008-05-28 


u9zs080 


3 


4 


m 


WFPC2/PC1/F555W 


400/700 


23.188 ±0.032 


23.205 ± 0.032 


2008-05-28 


u9zs080 


1 


2 


in 


WFPC2/PC1/F450W 


500/500 


22.994 ± 0.047 




2008-05-29 


u9zs090 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.109 ±0.029 


23.126 ±0.029 


2008-05-29 


u9zs090 


1 


2 


m 


WFPC2/PC1/F450W 


500/500 


22.836 ± 0.041 




2008-05-30 


u9zsl00 


1 


2 


in 


WFPC2/PC1/F450W 


500/500 


22.923 ± 0.048 




2008-05-30 


u9zsl00 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.075 ± 0.029 


23.092 ± 0.029 


2008-05-31 


u9zsll0 


1 


2 


in 


WFPC2/PC1/F450W 


500/500 


23.014 ±0.050 




2008-05-31 


u9zsll0 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.215 ±0.032 


23.232 ± 0.032 


2008-06-01 


u9zsl20 


3 


4 


m 


WFPC2/PC1/F555W 


400/700 


23.135 ±0.033 


23.152 ±0.033 


2008-06-01 


u9zsl20 


1 


2 


in 


WFPC2 /PC1/F450W 


500/500 


23.049 ± 0.046 




2008-06-02 


u9zsl30 


3 


4 


in 


WFPC2 /PC1/F555W 


400/700 


23.251 ±0.033 


23.268 ± 0.033 


2008-06-02 


u9zsl30 


1 


2 


in 


WFPC2/PC1/F450W 


500/500 


23.085 ± 0.048 




2008-06-03 


u9zsl40 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.132 ±0.038 


23.149 ±0.038 


2008-06-03 


u9zsl40 


1 


2 


in 


WFPC2/PC1/F450W 


500/500 


22.980 ± 0.043 




2008-06-04 


u9zsl50 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.261 ± 0.034 


23.278 ± 0.034 


2008-06-04 


u9zsl50 


1 


2 




WFPC2/PC1/F450W 


500/500 


22.977 ± 0.051 




2008-06-05 


u9zsl60 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.182 ±0.031 


23.199 ±0.031 


zUUo-UD-UO 


u9zsl60 


1 


2 


m 


VV r rVjZ [ r\^L/r VV 


oUU/ oUU 


zz.yUo it U.Uoz 




2008-06-06 


u9zsl70 


1 


2 


in 


WFPC2/PC1/F450W 


500/500 


23.073 ± 0.051 




2008-06-06 


u9zsl70 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.196 ±0.031 


23.213 ±0.031 


2008-06-07 


u9zsl80 


3 


4 


in 


WFPC2 /PC1/F555W 


400/700 


23.153 ±0.033 


23.170 ±0.033 


2008-06-07 


u9zsl80 


1 


2 


in 


WFPC2 /PC1/F450W 


500/500 


23.149 ± 0.053 




2008-06-08 


u9zsl90 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.247 ± 0.034 


23.264 ± 0.034 


2008-06-08 


u9zsl90 


1 


2 


m 


WFPC2/PC1/F450W 


500/500 


23.030 ± 0.056 




2008-06-09 


u9zs200 


3 


4 


in 


WFPC2/PC1/F555W 


400/700 


23.253 ±0.037 


23.270 ±0.037 


2008-06-09 


u9zs200 


1 


2 


m 


WFPC2/PC1/F450W 


500/500 


22.953 ± 0.043 




NGC 2403 X-l 2005-10-17 


j9h803020 




ACS/WFC /F606W 


1248 


24.564 ± 0.064 




2005-10-17 


j9h803010 




ACS/WFC/F435W 


1248 


24.687 ± 0.040 




2005-10-17 


j9h8a3010 




ACS/HRC /F330W 


2912 


23.243 ±0.171 





NGC 4559 X-7 



2001-05-25 
2001-05-25 
2001-05-25 
2005-03-08 
2005-03-08 
2005-03-08 
2005-03-09 



u6ch050 
u6ch050 
u6ch050 
j92w01040 
j92w01020 
j92w01010 
j92w02030 



9 


a 


b|c 


5 


6 


7|8 


1 


2 


3 4 



WFPC2/PC1/F814W 
WFPC2/PC1/F450W 
WFPC2/PC1/F555W 
ACS/HRC /F814W 
ACS/HRC /F555W 
ACS/HRC/F435W 
ACS/WFC /F606W 



500/500/500/500 
500/500/500/500 
500/500/500/500 

2400 

2400 

2400 

1828 



23.139 
22.878 
23.002 
22.981 
22.888 
22.744 
22.929 



±0.038 
±0.023 
±0.020 
± 0.083 
± 0.033 
±0.019 
±0.102 



23.016 ±0.020 
22.902 ± 0.033 
22.938 ±0.108 



NGC 5204 X-l 



2001-05-28 

2001- 05-28 

2002- 10-28 
2002-10-29 
2008-08-08 
2008-08-08 
2008-08-10 
2008-08-10 
2008-08-13 
2008-08-13 



u6723901r 
u6723902r 
j8dq01010 
j8dq01020 
u9zs210 
u9zs210 
u9zs220 
u9zs220 
u9zs230 
u9zs230 



1 


2 


in 


3 


4 


in 


1 


2 


in 


3 


4 


in 


1 


2 


in 


3 


4 


in 



WFPC2/WF2/F606W 600 22.492 ± 0.044 22.393 ± 0.089 

WFPC2/WF2/F814W 600 22.925 ±0.103 

ACS /HRC/F435W 2600 22.271 ±0.067 

ACS /HRC/F220W 2720 19.851 ± 0.024 

WFPC2/PC1/F450W 500/500 22.301 ± 0.022 

WFPC2/PC1/F555W 500/500 22.572 ± 0.021 22.601 ± 0.021 

WFPC2/PC1/F450W 500/500 22.350 ± 0.023 

WFPC2/PC1/F555W 500/500 22.528 ± 0.019 22.545 ± 0.019 

WFPC2/PC1/F450W 500/500 22.307 ± 0.023 

WFPC2/PC1/F555W 500/500 22.463 ± 0.018 22.478 ± 0.018 



NGC 5408 X-l 



2000-07-04 
2000-07-04 
2009-04-04 
2009-04-04 
2009-04-04 
2009-04-04 



u6724101r 

u6724102r 

ubah290[l|2|3]m 

ubah290[4|5|6|7]m 

ubah290[bm|dm|fn] 

ubah290 [cm | en | gn] 



WFPC2/WF3/F606W 
WFPC2 /WF3/F814W 
WFPC2/WF4/F439W 
WFPC2 /WF4 /F336W 
WFPC2 /WF4 /F555W 
WFPC2 /WF4 /F814W 



600 

600 
600/600/600 
1100/1100/1100/1100 
180/180/160 
180/180/160 



22.175 ±0.020 
22.266 ± 0.040 
21.879 ±0.028 
20.216 ±0.016 
22.114 ±0.022 
22.339 ± 0.045 



22.161 ±0.041 



22.132 ±0.022 



NGC 6946 ULX-1 



1996-01-27 
1996-01-27 
2001-06-08 
2001-06-08 
2001-06-08 
2004-07-29 



u33d0105t WFPC2/PC1/F555W 

u33d010[6 7]t WFPC2/PC1/F439W 

u6ch020 5 6|7|8]r WFPC2/PC1/F450W 

u6eh020[l 2|3|4]r WFPC2/PC1/F555W 

u6eh020 9 a|b|c]r WFPC2/PC1/F814W 

j8mxd8sxq ACS/WFC /F814W 



400 
400/400 
500/500/500/500 
500/500/500/500 
500/500/500/500 

120 



21.268 ±0.115 
21.239 ±0.184 
20.905 ±0.130 
21.238 ±0.122 
21.494 ±0.103 
21.450 ±0.072 



21.275 ±0.115 



21.267 ±0.122 



Note. — "ipiitcr is the magnitude quoted in the filter bandpass, mv is the magnitude in Johnson V , translated from F555W or F606W. All 
magnitudes are in VEGA zero point corrected for extinction. 



